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ABSTRACT 

£4T) , Correlations between the intrinsic shapes of galaxy pairs, and between the intrinsic 

If) ■ shapes of galaxies and the large-scale density field, may be induced by tidal fields. 

These correlations, which have been detected at low redshifts (z < 0.35) for bright 
red galaxies in the Sloan Digital Sky Survey (SDSS), and for which upper limits exist 
' for blue galaxies at z ~ 0.1, provide a window into galaxy formation and evolution, 

\ and are also an important contaminant for current and future weak lensing surveys. 

Measurements of these alignments at intermediate redshifts (z ~ 0.6) that are more 
relevant for cosmic shear observations are very important for understanding the origin 
k> , and redshift evolution of these alignments, and for minimising their impact on weak 

; i ■ lensing measurements. We present the first such intermediate-redshift measurement 

d ' for blue galaxies, using galaxy shape measurements from SDSS and spectroscopic red- 

shifts from the WiggleZ Dark Energy Survey. Our null detection allows us to place 
upper limits on the contamination of weak lensing measurements by blue galaxy in- 
trinsic alignments that, for the first time, do not require significant model-dependent 
extrapolation from the z ~ 0.1 SDSS observations. Also, combining the SDSS and 
WiggleZ constraints gives us a long redshift baseline with which to constrain intrinsic 
alignment models and contamination of the cosmic shear power spectrum. Assuming 
that the alignments can be explained by linear alignment with the smoothed local den- 
sity field, we find that a measurement of as in a blue-galaxy dominated, CFHTLS-like 
survey would be contaminated by at most io'03 (95 per cent confidence level, SDSS 
and WiggleZ) or ±0.03 (WiggleZ alone) due to intrinsic alignments. We also allow 
additional power-law redshift evolution of the intrinsic alignments, due to (for exam- 
ple) effects like interactions and mergers that are not included in the linear alignment 
model, and find that our constraints on cosmic shear contamination are not signifi- 
cantly weakened if the power-law index is less than ~ 2. The WiggleZ sample (unlike 
SDSS) has a long enough redshift baseline that the data can rule out the possibility 
of very strong additional evolution. 

Key words: cosmology: observations - gravitational lensing - large-scale structure 
of Universe - galaxies: evolution. 
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1 INTRODUCTION 

Gravitational lensing, the deflection of light due to matter 
between the source and the observer, is sensitive to all mat- 
ter (including dark matter) . As a result, in the past decade , 
weak gra v itatio nal lensing (|Bartelmann fc Schneider! l200ll ; 
iRefregieri l2003tl has become a powerful tool for ad- 
dressing outstanding questions related to cosmology and 
galaxy formation. Its scientific applications include mea- 
surement of the amplitude of matter fluctuations using 
the auto-correlation of galaxy shapes (e.g . , most recently, 
Hoekstra et al. 2006 ; Semboloni et al.l 2006; Benj amin et al.l 



20071 : iMassev et al.ll2007bl : ISchrabback et alJl200gh . known 
as cosmic shear; and determination of the relation- 
ship between the baryonic content and the dark matter 
content of galaxies, using the cross-correlation between 
background galaxy shapes and f o reground galaxy posi- 
tions (e.g. 



Hoekstra et~al] 120051 ; iHevmans et~aH l2006al ; 



iMandelbaum et al.ll2006ch . known as galaxy-galaxy lensing 
Because of the utility of these applications of lensing, plus its 
potential to constrain models of dark energy by splitting the 
sample o f source galaxi es into redshift slices (tomography: 
[Hull2002l ; lHutereril2002i ). future surveys are being planned 
to measure the lensing signal with sub-per cent statistical 
errors. 

There is a large body of work devoted to solving the 
technical problems in measuring the weak le nsing signal, pri- 
marily related to unbiased shear estimation (IHevmans et al.l 
l2006bl ; IMassev et alJl2007al: iBridle et al.ll2009f) and to pho- 



tometric redshifts dBernstein fc Jainl 20041; llshak fc Hiratal 



l2005l ; lHuterer et alj|2006l ; lAbdalla et alJl200Sl) . In this work, 
we focus on a source of astrophysical uncertainty, intrinsic 
alignments of galaxy shapes. When measuring the lensing 
signal, it is assumed that in the absence of lensing, galaxy 
shapes are uncorrelated. Intrinsic alignments are alignments 
of galaxy shapes that violate that assumption, for example 
due to the alignment of galaxy shapes with a local tidal 
field. These alignments can therefore contaminate the grav- 
itational lensing signal. 

One type of intrinsic alignment is the correlation be- 
tween the intrinsic ellipticities of two galaxies (II cor- 
relations) that reside in the same local or large-scale 
structure. Cosmological A-body simulations robustly pre- 
dict alignments between the shapes of dark matter halos 
that a re a declining func t ion of separation JSplinter et all 
19971; lOnuora fc Thomas! |2000|; iFaltenbacher et all l2002fc 



Hopkins et al.ll2005l ; ILee et al.ll200Sl ). However, the true ob- 
servational impact of these alignments are difficult to esti- 
mate using Af-body simulations, because the observed align- 
ments depend on the shape of the baryonic component of the 
galaxy rather than on dark matter alone. While several ana- 
lytical models for these alignments have also been developed 
((Croft fc Metzleri 120001; iHeayens et alJl200Cl : ICatelan et ail 
l200ll ; ICrittenden et al.ll200ll ; |jind[2002h . they predict wildly 
varying levels of alignment, so observational constraints are 
necessary. 

More recently, iHirata fc Seliakl (|2004 l pointed out that 
the correlation of galaxy shapes with large-scale density 
fields can also contaminate lensing measurements. These 
alignments, known as GI correlations, are caused by a lower 
redshift tidal field that both causes gravitational shear expe- 
rienced by a higher redshift galaxy, and intrinsically aligns 



the shape of galaxies that are in the tidal field. GI cor- 
relations are also predicted to have very different mag- 
nitudes depending on the model used to estimate them 
llHui fc Zhand l2002i ; IHirata fc Seliakl |2004| ; IHevmans et~aH 
l2006d\ and have been detected using dark matter ha- 
los at many different mass scales in A-body simulations 
jBailin fc Steinmet j|2005l;lAltav et al.ll2006l;lBasiIakos et all 
120061 ; IHevmans et al.ll2006d : iKuhlen et al.ll2007T l. 

The relevant signature of these GI correlation detec- 
tions in A-body simulations is that dark matter halos align 
so that they point preferentially towards other halos that 
are part of the same large-scale structure. When consider- 
ing GI correlations of galaxies that contaminate cosmolog- 
ical weak lensing measurements, the effect is manifested as 
galaxy shapes that point preferentially towards other galax- 
ies (both locally, within a halo, and on cosmological scales). 
These alignments of galaxy shapes are anti-correlated with 
the gravitational shear due to large-scale structure, so GI 
correlations reduce the measured cosmic shear signal, unlike 
the II correlations which increase it. Also unlike the II cor- 
relations, the GI correlations are not due to the inclusion of 
galaxy pairs at the same redshift; pairs at different redshifts 
are affected when the higher redshift galaxy of the pair is 
lensed by a structure that has caused an intrinsic alignment 
of the lower redshift galaxy. 

Several different schemes have been proposed to remove 
intrinsic alignment contamination from weak lensing mea- 
surements, including the removal of galaxy pairs that are 
close i n redshift space (to rem ove I LlKing fc Schneiddl2002l . 
120031 ; IHevmans fc Heavens! 120031 ; iTakada fc White! 120041 ); 



projecting out both types of intrinsic alignments using 
t heir known scalings with the redshifts of the galaxy pair 
llHirata fc Seliakl |2004|; Ijoachimi fc Schneider] |200S| ; IZhand 
120081 ; iJoachimi fc Schneider! |2009| )~ and modeling them 
jointly with the lensing signal using some parametric mod- 
els, the parameters of which are then marginalised over 
|Kinel2005l ; lBridle fc King|2007h . A common feature of these 
methods is a loss of information, and therefore weakening 
of cosmological constraints from the weak lensing signal. 
Measurements of intrinsic alignments can place strong pri- 
ors on the intrinsic alignment model, which would minimise 
the loss of cosmological information from future surveys. 
Direct intrinsic alignment measurements will also constrain 
the impact of intrinsic alignments on previous lensing mea- 
surements that did not explicitly account for them. This 
measurement is particularly important given the aforemen- 
tioned difficulty in theoretical predictions; however, the ob- 
servations can then be used to refine the theory and, in turn, 
learn something about galaxy formation and evolution. 

To observe intrinsic alignments, we require a source of 
data with robust galaxy shape measurements free of con- 
tamination from the PSF, and a way of isolating nearby 
(in all three dimensions) galaxy pairs. The GI correlations 
are then measured by calculating, statistically, the tendency 
for galaxies to point towards other galaxies that are rela- 
tively nearby (on tens of h~ x Mpc scales). Alternatively, it 
is possible to measure II correlations at low redshift without 
any redshift information, given that t he cosmic shear sig - 
nal below z ~ 0.2 is vanishingly small l|Brown et alj [20021 ). 
On the large scales used for cosmological lensing analyses, 
the first measurement of G I correlations used SDSS data 
jMandelbaum et all [2006b! ). with a follow-up analysis by 
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iHirata et ail ( |2007f ) that also included redshifts of SDSS 
Luminous Red Galaxies (LRG s) and from the 2d F-SDSS 
LRG and QSO survey C2SLAQ. ICannon et alJkood ) to con- 
strain intrinsic alignments of red galaxies up to intermediate 
redshifts, z ~ 0.5-0.6. While GI correlations were detected 
in these works at z ~ 0.1-0.4 for bright red galaxies (and 
II correlations for th e same galaxy sample were found by 
lOkumura et~al]|2009l ). with a weak (2<r) detection at inter- 
mediate redshifts (due to the small size of the 2SLAQ sam- 
ple), further work at intermediate to high redshift is cru- 
cial for constraining the impact of intrinsic alignments on 
cosmological lensing analyses. It is difficult to extrapolate 
these low-redshift analyses to higher redshift, because differ- 
ent dynamical scenarios might entail very different redshift 
evolution. For example, blue galaxies at z ~ 0.1 that have 
no measurable GI alignment in SDSS may have been very 
highly aligned with the density field at intermediate-high 
redshift, with mergers and interactions serving to disrupt 
those alignments, leading to the null detection that we see 
at low redshift; or, these alignments might be very small at 
all redshifts. 

Due to its overlap with the SDSS, which provides 
gala xy shape measurements, the WiggleZ Dark Energy sur- 
vey (|Drinkwater et al"1l2010h is an ideal source of spectro- 
scopic redshifts for constraining intrinsic alignments of UV- 
selected blue galaxies at intermediate redshift. In this work, 
we use that sample to attempt the first measurement of 
galaxy intrinsic alignments for blue galaxies at intermediate 
redshift, which fills in a very important gap in our knowl- 
edge of intrinsic alignments. At higher redshift, we expect 
that blue galaxies will dominate the galaxy samples used 
for weak lensing. Thus, these observations will facilitate fur- 
ther development in the fields of weak lensing and galaxy 
dynamics and evolution. 

Here we note the cosmological model and units used 
throughout this work. Pair separations are measured in co- 
moving ft -1 Mpc (where Ho = 100/ikms -1 Mpc" 1 ), with 
the angular diameter distance computed in a spatially flat 
ACDM cosmology with Q m = 0.3. For the bias and cos- 
mic shear calculations, we additionally normalise the mat- 
ter power spectrum using erg = 0.75, set the baryon density 
fib = 0.05 and scalar primordi al spectral index n 3 = 1, and 
use the transfer function from iMal (|l996h . 

We begin in Section [2] with a summary of the intrinsic 
alignment and cosmic shear formalism used in this work. 
Section |3] contains descriptions of data used for the analysis. 
The methodology used for the data analysis is described in 
Section 21 We present the results of the analysis in Section 
including systematics tests and a comparison with previous 
observations. The interpretation of these results, including 
an estimate of contamination of the cosmic shear signal, is 
given in Section |SJ and we conclude in Section [7] 



ponents: the gravitational lensing- induced shear <y , and the 
"intrinsic shear" -y 7 , which includes any non-lensing shear, 
typically due to local tidal fields. Therefore we can write 
the U-mode shear power spectrum between any two red- 
shift bins a and j3 as the sum of the gravitational lens- 
ing power spectrum (GG), the intrinsic-intrinsic, and the 
gravitational-intrinsic terms, 



C EE (a/3) = cf E ' GG 



(ap) + C BEJ1 (aP) + C? E ' GI (aP). (1) 



iMandelbaum et al.l (|2006bl ) presented the Limber integrals 
that allow us to determine each of these quantities in terms 
of the matter power spectrum and intrinsic alignments 
power spectrum. In a flat universe, the GI contamination 
term can be written as 



^yEE.GI f o\ 

Cj («P) 



-3U(r)Wp(r)Ps* t [k = 



1 + 1/2 



(2) 



where rn is the comoving distance to the horizon, f a (r) is 
the comoving distance distribution of the galaxies in sample 
a, and 



W a (r) = ^ft m H^(l + z) 



r(r' 



-U(r')dr'. 



(3) 



The generalisat ion of these equations to curved universes 
can be found in IMandelbaum et al.ll2006bl . 

The density-intrinsic shear cross-power spectrum 
P s (k) that enters into Eq. ((2J is defined as follows. If one 
chooses any two points in the SDSS survey, their separation 
in redshift space can then be identified by the transverse 
separation r v and the radial redshift space separation IT. 
The + and x components of the shear are measured with 
respect to the axis connecting the two galaxies (i.e., positive 
+ shear is radial, whereas negative + shear is tangential). 
Then one can write the density-intrinsic shear correlation in 
Fourier space as 



(4) 



P St¥ {k) = -2n / & + (r p ,n)J 2 (kr p )r p dr p dn, 



where £a+(r" p , IT) is the correlation function between the den- 
sity contrast S = p m /~p m — 1 and the galaxy density- weighted 
intrinsic shear, 7+ = (1 + <5 S )7 J (where S g = p g /~p g — 1). It 
is often convenient to do the projection along the radial di- 
rection, 



w s+ (r p )= / c«+(r P ,n)dn. 



(5) 



A similar set of equations can be written for the 
intrinsic-intrinsic terms. For example, we can define 

Cf E ^(a,3) = r %Ur)U{r)P?F{k = l -±l£, r) (6) 



2 FORMALISM 

Here we briefly summarise the formalism for the analysis 
of intrinsic alignment contamination to the lensing shear 
co rrelation function. Ou r notation is consistent with that 
of IHirata fc Seliakl (|2004f ), IMandelbaum etafl (|2006bh . and 



IHirata et all <|2007l ) 



The observed shear 7 of a galaxy is a sum of two com- 



in terms of the E-mode power spectrum of the density- 
weighted intrinsic shear, P EE ■ Likewise, the intrinsic- 
intrinsic correlations are 

P°i E (k) = f[^+ + (r p ,U)J+(kr p )+^ x (r p ,U)J-(kr p )] 

X 27rrp dr p dn (7) 

where J±(x) = [Jo(x) ± Ja{x)]/2. 

While (to first order) the lensing shear does not induce 
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any B-mode signal, intrinsic alignments are one of several ef- 
fects t hat may contribute to a nonzero _B-mode pow er spec- 
trum (|Hirata fc Seliakl |2004 iHevmans et al.1 l2006cT l. How- 
ever, we will not focus on these effects here. 



3 DATA 
3.1 WiggleZ 

The Wiggl eZ Dark Energy Survey at the Anglo- Australian 
Telescope (|Drinkwater et aTl I2010T 1 is a large-scale galaxy 
redshift survey of bright emission-line galaxies mapping a 
volume of order 1 Gpc 3 over the redshift range z < 1. The 
survey, which began in August 2006 and is scheduled to 
finish in July 2010, is obtaining ~ 200,000 redshifts for UV- 
selected galaxies covering ~ 1000 deg 2 of equatorial sky. It 
is performed using the multi-fibre spectrograph AAOmega, 
which can simultaneously obtain spectra for up to 392 
galax ies over a 2-degree-diameter field-of-view (|Sharp et al.l 
I2006T I. The principal scientific goal is to measure the baryon 
acoustic oscillation signature in the galaxy power spectrum 
at a significantly higher redshift than existing surveys. The 
target galaxy population is selected from UV imaging by 
the Galaxy Evolution Explorer (GALEX) satellite, matched 
with optical data from the Sloan Digital Sky Survey (SDSS) 
and Red Cluster Sequence survey (RCS2) to provide accu- 
rate positions for fibre spectroscopy. 

In this paper, we analyse the subset of the WiggleZ sam- 
ple lying in the SDSS survey areas assembled up to the end 
of the 09 A semester (May 2009). Specifically, we include 
data from the WiggleZ 9-hr (09h), 11-hr (llh), and 15-hr 
(15h) regions, centred at the following positions: (141.3, 3.6), 
(162.5, 3.5), and (220.0, 2.0) respectively (all positions are in 
degrees of right ascension and declination, J2000 equatorial 
coordinates). These regions together include 76 084 galax- 
ies with spectroscopic r edshifts classified as rel iable (with 
quality Q = 3, 4, 5, see iDrinkwater et alJ [2oToh . with the 
three regions containing 22 011, 21 746, and 32 327 galax- 
ies, respectively. The number density of galaxies with suc- 
cessful redshift estimates varies within these regions from 
~ 200 to ~ 340 per square degree, depending on how com- 
pletely a given region was observed. These galaxies consti- 
tute an e xtended sample rel ative to the one used for the 
analysis in Blake et al. (2009). As a consequence of the con- 
tinuing GALEX imaging campaign in the WiggleZ survey 
regions, the UV magnitudes of a fraction of the targets 
have been refined as the survey progresses, causing some 
originally-observed galaxies to now fail the survey magni- 
tude and colour selection cuts. This subset of galaxies was 
not in cluded in the original clustering analysis of lBlake et al.l 
l|2009l ). but has now been accommodated following suitable 
modifications to the random catalogue generation proce- 
dure. Full details will be presented in a future paper. 

The redshift error rate for galaxies that are assigned a 
redshift, whic h is a functi o n of r edshift, is given in Table [T] 
As shown in iBlake et al.1 J200I) figure 13, the typical B- 
band luminosity of this sample ranges from two magnitudes 
below Lj, at z = 0.2, to L« at 2 = 0.5, to two magnitudes 
above L, at z = 1. For simpler comparison against SDSS, 
Table [T] shows the rest-frame r-band magnitudes for this 
sample. These were derived on average from the UV and r 



band photometry, using a Lyman Break Galaxy template 
that is consistent with the WiggleZ galaxies being detected 
in the near-UV (NUV) but being far-UV dropouts. This 
template is a constant star formation rate model, with sig- 
nificant dust extinction added to match the ob served NUV- 
r colo ur versus redshift relation. Figure 3 of IWvder et al.l 
|2007l ). which shows the NUV-r colour-magnitude relation 
for galaxies with measurements from GALEX and SDSS 
around z = 0.1, is a good illustration of the nature of the 
WiggleZ sample. That figure shows a distinct, well-defined 
red sequence and a blue cloud, where the WiggleZ selection 
of NUV-r < 2 picks out the very blue edge of the blue cloud. 

We create random catalo gues for each Wig gleZ region 
using the method described bv lBlake et al.l (|2009f ) with mod- 
ifications to account for the use of an extended sample. In 
brief, random realisations of "parent" catalogues are first 
created which trace the variation in WiggleZ target density 
with Galactic dust extinction and GALEX exposure time. 
These realisations are then processed into random "redshift" 
catalogues by imposing the observing sequence of telescope 
pointings. The fraction of successful redshifts in each point- 
ing varies considerably depending on weather conditions. 
Furthermore, the redshift completeness within each point- 
ing exhibits a significant radial variation due to acquisition 
errors at the plate edges, which is also modelled. 

In order to model the observed intrinsic alignment sig- 
nal, we require a measurement of the bias of the sample used 
to trace the density field, which is the full WiggleZ redshift 
sample (including those galaxies without shape measure- 
ments). The galaxy biases that we use for thi s analy sis were 
measured using the method of IBlake et alJ (|2009h . which 
(in brief) involves the following steps: (i) measurement of 
the correlation function projected along the line of sight to 
±20/i _ Mpc; (ii) fitting this measurement to a power-law 
correlation function (also integrated along the same line- 
of-sight range, and with a model for redshift space distor- 
tions) to determine a correlation length; (iii) generating a 
dark matter correlation function using CAMB (|Lewis et al.l 
120001 ) with halofit = 1, assuming Q. m = 0.3, Q.b/0, m — 0.15, 
h = 0.7, n s = 1, and erg = 0.9, to determine the dark mat- 
ter correlation length; and (iv) estimating a galaxy bias us- 
ing the ratio of the correlation lengths, while accounting for 
the linear growth factor. Since we use as = 0.75 for the 
cosmic shear po wer spectrum calcul ations in this paper, we 
also increase the lBlake et al.l (|2009l ) bias measurements by 
(0.9/0.75) so that they are consistent with this cosmological 
parameter choice. Finally, we correct for the effect of red- 
shift blunders by decreasing theoretical model predictions 
by a factor of (1 — /bad) 2 , which corresponds to increasing 
the measured correlations by 1/(1 — /bad) 2 - 

The observed redshift distribution is, to good approxi- 
mation, a double Gaussian, with 77.6 per cent of the galaxies 
in a Gaussian with mean (z) = 0.595 and width a z = 0.236, 
and the remaining 22.4 per cent in a narrower Gaussian 
with mean (2) = 0.558 and width a z = 0.112. We will use 
this description of the redshift distribution in our theoretical 
modeling of the observations. 

3.2 SDSS spectroscopic sample 

In this paper, we compare the WiggleZ results with previous 
measurements of intrinsic alignments using az ~ 0.1 sample 
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Table 1. Parameters of the WiggleZ and SDSS samples used here: minimum redshift, maximum redshift, pair-weighted effective rcdshift, 
bias of the density field sample corrected for the fraction of bad redshifts, average absolute magnitude of the shape-selected and density 
field samples (averaged over 10°' 4M ) A;-corrected to the rest-frame SDSS r-band, fraction of bad redshifts, and fraction of WiggleZ 
galaxies with good shape measurements. Values used for the comparison with SDSS low-rcdshift blue galaxies are also provided for 
completeness. The magnitude values are M — 51og/i. 



Sample 




^max 


Zeff 


^density 


-^shape 


-^density 


/bad 


,/match 


WiggleZ all 


0.01 


1.3 


0.51 


1.50 ±0.04 


-20.9 


-20.7 


0.050 


0.33 


WiggleZ z < 0.52 


0.01 


0.52 


0.37 


1.28 ± 0.04 


-19.9 


-19.4 


0.032 


0.34 


WiggleZ z > 0.52 


0.52 


1.3 


0.62 


1.63 ± 0.04 


-21.2 


-21.0 


0.064 


0.32 


SDSS Blue L4 


0.02 


0.19 


0.09 


1.12 ±0.04 


-20.8 


-20.8 


< 0.01 


0.90 



of SDSS spectroscopic galaxies that are blue a nd have lumi- 
nositie s near L„ , the "blue L4" sample from iHirata et al.l 
(|2007t ). The properties of this sample are also given in Ta- 
ble [T] In that paper, we were interested in robustly iso- 
lating the red sequence from the blue cloud, so the colour 
separator used there defined this "blue" sample to include 
the entire blue cloud. The r-band luminosities used to de- 
fine this sa mple were fc-correcte d to z = 0.1 using KCOR- 
rect v3_2 (jBlanton et al.l I2003T ) with Petrosian apparent 
m agnitudes, extinction corrected using the reddening maps 
of ISchlegel et al.l (1 1998 ) and extinction-to-reddening ratios 
from IStoughton et al.l ( 20021 '). For fair comparison with the 
WiggleZ sample, here we calculate luminosities with the 
model magnitudes, and correct to z = 0. 



galaxies must (a) have extinction-corrected r-band model 
magnitude m r < 21.8, and (b) be well-resolved compared 
to the PSF size in both r and i bands (as quantified by 
the adaptive moments of the PSF and galaxy image). The 
WiggleZ galaxies with shape measureme nts were part of the 
genera l SDSS shape catalog presented in lMandelbaum et al.l 
(|2005t l and used for many subsequent science papers. Thus, 
they have already been subjected to all systematics tests 
detailed in those papers, p articularly the original paper and 
iMandelbaum et alj ()2006ah . which has other significant tests 
of shear systematics. Nonetheless, in this paper we will still 
present additional systematics tests to rule out the possi- 
bility that this sample has some unusual set of systematics 
compared to the rest of the shape catalog. 



3.3 Galaxy shape measurements 

For a subset of the WiggleZ galaxies, w e use shape mea - 
surements from the SDSS. The SDSS (|York et all |2000| ) 
imaged roughly n steradians of the sky, and followed 
up approximately one million of the d etected objects 
spect r oscopicallv dEisenstein et al ] 1200 ll : iRichards et al] 
l2002i ; IStrauss et al] I2002T ). The imaging was carried 
out b y drift-scanning the sky in phot ometric condi- 
tions Irlogg et al] l200ll: llvezic et al] l2004h in fi ve bands 



using 



(ugriz) (|Fukugita et al] Il99rj ; ISmith et alj 12002!') 
specially-designed wide-field camera (|Gunn et ail Il998f ). 
These imaging data were used to create the galaxy shape 
measurements that we use in this paper. All of the 
data were processed by completely automated pipelines 
that detect and measure photometric pro perties of ob- 
jects , and astrometric a lly calibrate the da ta (|Lupton et al] 
l200ll ; IPier et al.l 120031 ; iTucker et al1l2006l ). The SDSS has 
had seven major data releases, and is now complete 
dStoughton et al 1 12002|; lAbazaiian et al] 120031. 12004 1 2005} 
Finkbeiner et al] 12004: I Adelman- McCarthy et al] boOfj . 



20071 . l200Sl : lAbazaiian et alj|2009r i. 

We use the galax y ellipticity measurements by 
IMandelbaum et al] (|2005t ). who obtained shapes for more 
than 30 million galaxies in the SDSS imaging data down to 
extinction-corrected magnitude r = 21.8 usi ng the Reglens 
pipeli ne. We refer the interested reader to IHirata fc Seliakl 
l|2003h for an outline of the PSF correction t echni que (re- 
Gaussianization) and to IMandelbaum et al] (120051 ) for all 
details of the shape measurement. The full details of re- 
strictions imposed on ga laxy shape measurements are in 
IMandelbaum et al] (|2005t ) , but the two main criteria for the 
shape measurement to be considered high quality are that 



In the 09h, llh, and 15h fields, the fractions of WiggleZ 
galaxies with high-quality shape measurements are 33, 34, 
and 32 per cent, respectively, giving a number density that 
ranges from 60 to 100 degree" 2 . Fi gure[T] shows the redshift 
and r-band model magnitude distributions for the WiggleZ 
redshift sample, and the fraction with good shape measure- 
ments. As shown, the probability that there is a good shape 
measurement exhibits a significant magnitude dependence, 
and goes to zero for r > 21.8 due to a cut imposed on 
the shape catalogue. However, the redshift distribution of 
those galaxies with shapes is not substantially different from 
that of the full sample, because at any given redshift the 
most luminous galaxies tend to be large enough in apparent 
size relative to the PSF that they have a measurable shape. 
There is also a slight region-to-region variation of the match 
fraction, for two main reasons: because the typical seeing 
in the SDSS observations varies with position, and because 
the pointing strategy of the WiggleZ observations prioritises 
fainter galaxies, which are less likely to have a good shape 
measurement, so the more completely observed regions will 
tend to have a higher match fraction. 

To do this measurement, we also need random cata- 
logues that correspond to the shape-selected subset of the 
galaxies. To flag a fraction of galaxies in our random galaxies 
as possessing "good shapes," we estimated (from the data) 
the probability of a galaxy possessing a good shape as a 
function of the seeing of the SDSS observation and of the 
galaxy magnitude, and imposed this probability function on 
the random points (separately in each region). As shown in 
Fig. [T] the good shape fraction is a decreasing function of 
magnitude; as expected, we find that it is also a decreasing 
function of the seeing FWHM. 
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Figure 1. Top: Solid black, dashed red, and long-dashed blue 
lines show the fraction of WiggleZ galaxies in the 09h, llh, and 
15h fields (respectively) that have high-quality shape measure- 
ments in SDSS, as a function of r-band apparent model magni- 
tude. The arbitrarily-normalised, hatched magenta curve drawn 
with dotted lines shows the apparent magnitude distribution of 
the full WiggleZ sample. Bottom: Similar to the top, but as a 
function of redshift. The local minimum in the good shape frac- 
tion at z as 0.3 is created by the strong correlation of galaxy 
luminosity with redshift in the WiggleZ sample, with the result 
that galaxies at lower redshifts have preferentially smaller effec- 
tive radii, which more than offsets the larger apparent size due to 
the lower redshift. The tail of good shapes at z > 1 is partially 
due to the redshift blunder rate in the WiggleZ sample, which is 
about 3% at z = 0.6 but rises steeply at z > 1 to almost 50% 
iBlake et al.ll2009l . fig. 6). 



4 METHODOLOGY 

The software for comput ation of correlation functio ns is 
the same as that used in iMandelbaum et all (|2006bJ ) and 
iHirata et all |2007l) . In order to find pairs of galaxies, this 
code uses the SDSSpix packageQ To reduce noise in the 
determination of galaxy-random pairs, we use 100 random 
points for each real galaxy in the catalogue. The correla- 
tion functions are computed over a 120 /i _1 Mpc (comov- 
ing) range along the line of sight from II = —60 to II = 
+60h~ 1 Mpc, divided into 24 bins with size All = 5/i _1 Mpc, 
and the projected correlation function is computed by "inte- 
gration" (technically summation of the correlation function 
multiplied by All) over IT. This value of n max was chosen to 
minimise the loss of cor related galaxy pairs at all projected 
separations used here |Padmanabhan et al.l [20071 ) without 
increasing the noise excessively. We also show results with 
II max = 20/i _1 Mpc, which have better S/N but are more 
complicated to interpret due to redshift space distortions, 
and therefore are not used for cosmological interpretation in 
this paper. 



1 http : //lahmu.phyast .pitt . edu/~scranton/SDSSPix/ 



This calculation is done in iVbin = 10 radial bins from 
0.3 < r p < 60/i _1 Mpc. Covariance matrices are determined 
using a jackknife with 49 regions, in order to account prop- 
erly for shape noise, shape measurement errors, and cosmic 
variance. This number was chosen to be large enough to ob- 
tain a stable covariance matrix for the fits (it must be larger 
than see Appendix D of lHirata et al.|[2004l ) but small 

enough that the size of a given jackknife region is larger than 
the scale on which the correlation is to be measured. Each 
such calculation is carried out separately for the three Wig- 
gleZ regions, to check for consistency, before averaging over 
the regions. 

The code measures several different correlation func- 
tions simultaneously; here we describe the estimator for each 
one. 

For the GI cross-corre lation function £q+(r p , If), we use 
a generalisation of the LS (|Landv fc Szalav Il993l ) estimator 
for the galaxy correlation function. This generalisation can 
be expressed as 

S+(D-R) _ S+D-S+R 



£«+(»!>, n) = 



(8) 



R s R Rs R 

where S+D is the sum over all real ("data") galaxy pairs 
(with one galaxy in the subset with shapes, and the other 
galaxy in the full WiggleZ redshift sample) with separations 
r p and II of the + component of shear: 



2TZ : 



(9) 



where e+(j\i) is the + component of the ellipticity of shape 
sample galaxy j measured relative to the direction to den- 
sity field galaxy i, and TZ is the shear responsivity (that rep- 
resent s the response of our ellipticity definition to a small 
shear; iKaiser et all 1 19951 . iBernstein fc Jarvidliool ). S-\-R is 
defined by a similar equation, but using pairs derived from 
the real sample with shape measurements and the full ran- 
dom catalogues. R a R is the number of pairs of random galax- 
ies with separations r v and II such that one of those ran- 
dom galaxies is in the subset that is statistically likely to 
have a good shape measurement in SDSS, and the other is 
in the full WiggleZ random sample. (S+R and R S R are un- 
derstood to be rescaled appropriately since the number of 
random catalogue galaxies differs from the number of data 
galaxies.) Note that when doing the summation in Eq. ((9]l 
to determine S+D (or the comparable summations for S+R 
and R S R), we use all pairs regardless of which galaxy, i or 
j (in the density field tracer and intrinsic shear tracer sam- 
ples, respectively), is in the foreground. The reason for this 
choice is that we are attempting to detect an alignment due 
to the two galaxies experiencing the same tidal field since 
they are in close 3D proximity, rather than a lensing effect 
(which would require them to be at different redshifts). 

Averaged over a statistical ensemble, (S+) — (D — R) — 
0, so that systematics in the shear or the number density 
cancel to first order. Positive £ 9 + indicates a tendency to 
point towards over-densities of galaxies (i.e., radial align- 
ment, the opposite of the convention in galaxy-galaxy lens- 
ing that positive shear indicates tangential alignment). For 
the purpose of systematics tests, we can define an analogous 
estimator using the other ("x") ellipticity component for 
Csx ■ 

For the intrinsic shear auto-correlation functions 
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£ ++ (r p , IT) and £ x x (r P , II), we restrict ourselves to the sub- 
set ol the data with shape measurements, and use the esti- 
mators 



- S+S+ ~ S X S X 

C++ = o n and ?xx = 



where 
S4-S4- 



R s Rs 



E 



R s Rs 



e+(j\i)e+(i\j) 
(2K) 2 



(10) 



(11) 



(with both i and j denoting galaxies in the shape-selected 
sample) and similarly for S X S X . Since (S+) = (S x ) = 0, 
the cancellation of systematics to first order works again, i.e. 
the square of any spurious source of shear adds to Eq. (|10[1 
instead of the shear itself. Projected quantities such as w g + 
and w++ are then obtained by line-of-sight integration. 

Note that the intrinsic shear autocorrelation, w++, may 
potentially have some contribution from cosmic shear. How- 
ever, for the median redshift of the WiggleZ sample, the 
predicted contribution from cosmic shear to w++ (for a con- 
cordance co smology) is of ord er 10 -3 at 10 h~ x Mpc (see, 
for example. I Jarvis et al. I l2006l ; to estimate the cosmic shear 
contribution to u>++ we must include a factor of 2II max due 
to the line-of-sight integration). We shall see in Section [5] 
that this cosmic shear contribution is well below our error- 
bars and therefore undetectable. This is a consequence of 
the very low galaxy number density for the subset of the 
WiggleZ sample that has good shape measurements, which 
means that a cosmic shear measurement is not feasible (even 
were we to avoid the restriction that the pairs be close along 
the line-of-sight, which is necessary to detect intrinsic align- 
ments, but not cosmic shear). 



5 RESULTS 

We begin by presenting the projected intrinsic alignment 
cross-correlation functions, w g + and w++. The results are 
shown for each region in Fig. [21 We have scaled the signal 
by r p ' 8 for easy viewing. As shown, both w g + and are 
consistent with zero in all regions. Furthermore, there is no 
sign of any systematic discrepancy between the results in 
different regions, so for all subsequent tests we present only 
the results averaged over region (shown in Fig. Adja- 
cent radial bins on large scales are correlated at the level 
of several tens of per cent for w g +, whereas w++ is suffi- 
ciently dominated by shape noise that the points are nearly 
uncorrelated. 

As shown in Fig. [2] the errorbars for w g + and w++ are 
within a factor of two of each other. It is worth considering 
at this point what value the w++ measure ment has for con- 
strain ing intrinsic alignments. As shown in lHirata fc Seliakl 
(|2004 ). the linear alignment model at z ~ 0.5 predicts that 
for typical scales used in this measurement, the ratio of the 
II to the GI power spectra is of order 0.2. Thus, within the 
context of the linear alignment model, our non-detection of 
GI alignments implies an even lower II amplitude that, given 
our errorbars, will be undetectable. We can conclude that 
the non-detection of II does not give us significant additional 
information about intrinsic alignments within the context of 
the linear alignment model. The utility of our II measure- 
ment is that it allows us to (a) rule out significant shear 




1 10 
r p [fr'Mpe] 

Figure 2. Top: Projected GI cross-correlation signal w g +(r p ), 
multiplied by r^- 8 . Results are shown for each field separately: 09h 
(black solid line with hexagonal points); llh (red dashed line with 
crosses); and 15h (blue dot-short dashed line with solid squares). 
Points at a given value of r p are slightly horizontally offset for 
clarity. Bottom: Same as the top, but for the II cross-correlation 
signal w ++ (r p ). 



systematics that would lead to galaxy shape correlations 
significant enough that they might affect w 9 +, and (b) rule 
out substantial intrinsic alignments due to other causes be- 
sides the linear alignment model. For e xample, the simplest 
form of the quadra t ic alignment model ([Catelan et alj|200ll ; 
iHui fc Zhan3l2002l ; lHirata fc Seliakll2004l ) predictTzero GI- 
type alignments, but nonzero II alignments. 

To assess the consistency of these signals with zero, we 
include Table [2] This table shows the \ 2 f° r a fit 1° zero s ig~ 
nal, including correlations between radial bins (by using the 
full inverse covariance matrix). We also include the prob- 
ability for a random vector with this covariance matrix to 
exceed the given \ 2 by chance, p(> x 2 ). To calculate this 
probability value, we included the fact that the jackknife 
covariance matrices lead to a x 2 value that does not follow 
the expected distribution, because of noise from the finite 
number of jackknife regions. To include the effects of noise, 
we use a simulation ba sed on the formalism in Appendix 
D of Irlirata et alj (|2004 ). As shown in the first two lines of 
Table [2j the w g + and w++ shown in Fig. [2] are indeed con- 
sistent with zero, given that our criterion for inconsistency 
with zero is p(> \ 2 ) < 0.05. 

Next, we split the galaxies at z — 0.52 and recompute 
these correlation functions for each redshift subsample, with 
effective redshifts of 0.38 and 0.63, respectively (see Table[l|. 
This value of redshift was chosen to give approximately 1/3 
of the galaxies at z < 0.52, and 2/3 above, which (given the 
higher measurement noise in the latter sample) yields ap- 
proximately equal S/N for the w g + and w++ measur e ments 
in the two redshift slices. As shown in iBlake et all ([2009) 
and our Table 1, this split corresponds to a luminosity split. 
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Table 2. Comparison between our measurements and zero, using 
all ten radial bins, for the signals averaged over region. 



IT range 
(h- 1 Mpc) 



z range 



statistic 



X 2 P(> X 2 ) 





n| < 60 


All 


Wg+ 


4.83 


0.90 




n| sg 60 


All 


W++ 


14.53 


0.25 




n «; 60 


z ^ 0.52 


w g+ 


8.98 


0.60 




n ^ 60 


z ^ 0.52 


w ++ 


13.28 


0.30 




n < 60 


z > 0.52 


w 9 + 


7.24 


0.84 




n < 60 


z > 0.52 


W++ 


18.08 


0.13 




n ^ 20 


All 


w g+ 


15.82 


0.20 




n ^ 20 


All 


W++ 


14.92 


0.23 




n < 60 


All 


Wgx 


3.67 


0.96 




n < 60 


All 


1V+X 


18.21 


0.13 




n sg 60 


All 


»XX 


8.60 


0.62 


100 «: |n| =g 150 


All 


w 9 + 


8.70 


0.62 


ioo < n| s: 150 


All 


W + + 


5.60 


0.85 
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Figure 3. Top: Projected GI cross-correlation signal iu s _|_(r p ), 
multiplied by r p s . Results are shown averaged over all regions, 
for the two redshift subsamplcs. Points at a given value of r p are 
slightly horizontally offset for clarity. Bottom: Same as the top, 
but for the II cross-correlation signal w++(r p ). 



Fig. [3] shows that the results for the two redshift slices are 
consistent with each other, with no detection of any intrinsic 
alignment signal. These null results can also be confirmed 
using the third through sixth lines of Table [2] 

Finally, in Fig. [4] we show the results with the two differ- 
ent values of II max , 20 and 6O/1 -1 Mpc. As previously noted, 
the former results (while less noisy) are more complex to 
interpret due to the need for a model of redshift-space dis- 
tortions, and so we use the latter for all cosmological inter- 
pretation. However, we can confirm again in lines 7 and 8 
of Table [2] that the results with the smaller II max value are 
also consistent with zero. 



1 - 



-1 - 




_n_ = 60h~'Mpc 
- n m»x = 20h-'Mpc 




r p [fr'Mpe] 

Figure 4. Top: Projected GI cross-correlation signal w g +(r p ), 
multiplied by . Results are shown averaged over all regions, 
for the two values of II max . Points at a given value of r p are 
slightly horizontally offset for clarity. Bottom: Same as the top, 
but for the II cross-correlation signal iu_|_+(r p ). 



5.1 Systematics tests 

In this section we present several systematics tests, though 
the galaxy shape measurem ents used for this paper were 
already tested extensively in lMandelbaum et al.l ()2005l ) and 
subsequent papers. While these tests may seem irrelevant 
given the null results for w g + and we would like to 

rule out the possibility that a real astrophysical signal may 
be masked by a systematic error of similar magnitude but 
opposite sign. We also, however, note that the most likely 
sign of contributions of both PSF systematics and intrinsic 
alignments to w++ is positive, so the null detection itself 
constitutes a constraint on systematics. 

The first test involves the use of the other ellipticity 
component to compute w g x and w+ x (which should be zero 
by symmetry for a real astrophysical signal, since intrinsic 
alignments only induce alignments in the radial/tangential 
direction, but may be generated due to certain errors in PSF 
correction) . As shown in Fig. [5] there is no sign of either of 
these signals; they are completely consistent with zero, as 
confirmed in Table [2] 

The second systematics test is to compute these signals 
using pairs at large line of sight separations, 100 ^ |LT| < 
150ft -1 Mpc. This test will help show whether there is any 
spurious signal due to some systematic effect, since those 
pairs are effectively not correlated. Figure [6] shows no sign 
of nonzero w g + or uj++ for pairs at large line-of-sight sep- 
arations, indicating that the potential contaminants to the 
signal are within the errors. Table [2] confirms this finding 
quantitatively. 
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Figure 5. Top: Projected GI cross-correlation signal «J 9 + (r p ) and 
the systematics test w gx (r p ) (as indicated on the plot), multiplied 
by rp 1 - 8 . Results are shown averaged over all regions. Points at a 
given value of r p are slightly horizontally offset for clarity. Bottom: 
Same as the top, but for the II auto-correlation signals ui_|__|_(r p ) 
and w xx (r p ), and their systematics test, the cross-correlation 
w +x (r p ). 
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Figure 6. Top: Projected GI cross-correlation signal ui 3 _|_(r p ) 
for associated and non-associated galaxy pairs (as indicated on 
the plot), multiplied by rp 1 - 8 . Results are shown averaged over 
all regions. Points at a given value of r p are slightly horizontally 
offset for clarity. Bottom: Same as the top, but for the II cross- 
correlation signal w++(r p ). 



Figure 7. Projected GI cross-correlation signal w g +(r p ) multi- 
plied by rp 1 - 8 , for several galaxy samples as labelled on the plot. 
Results for WiggleZ are shown averaged over all regions. Points 
at a given value of r p are slightly horizontally offset for clarity. 



5.2 Comparison with previous observations 

There have been several previous measurements of large- 
scale intrinsic alignments. We focus on those that are pre- 
sented using comparable estimators that include the ellip- 
ticity (i.e., not those that correlate the position angles) and 
on those that go to the large scales that are of interest for 
cosmic shear. 

First, iMandelbaum et al.1 (|2006bh presented GI and II 
correlations for SDSS Main spectroscopic sample galaxies 
(typical z ~ 0.1) split into luminosity bins, which included 
a pos itive detection of G I signal for the bins with L > L,. 
Next. lrlirata et all (|2007l ) showed results for the Main sam- 
ple split into both colour and luminosity bins, in addition to 
new results from the SDSS LRG sample (red galaxies with 
z ~ 0.3) and from the 2SLAQ redshift survey using shape 
measurements from SDSS. 

On Fig. [7] we show the measured w g + for the Wig- 
gleZ sample, for one of the four blue galaxy samples derived 
from the SDSS Main galaxy sample (L4, with ( L) ~ L«); for 
the in termediate SDSS LRG luminosity bin in lHirata et alj 
(|2007l ); and for the 2SLAQ sample (a red galaxy sample a t 
the same typical redshift as WiggleZ, ICannon et al.1 120061 ). 
As shown in this figure, the signals for "typical" blue galaxies 
at z ~ 0.1 (SDSS blue L4), and for UV-selected galaxies at 
z ~ 0.6 (WiggleZ) are consistent with zero. The constraints 
are tighter using SDSS than using WiggleZ despite the larger 
volume of the WiggleZ survey, because the WiggleZ survey is 
shot-noise limited on small scales given the sparse sampling 
of the galaxies that are targeted for spectroscopy. However, 
the WiggleZ measurement is the first at cosmologically rel- 
evant redshifts for blue galaxies, which tend to dominate 
cosmic shear samples. 

Detailed comparison of the SDSS blue L4 sample and 
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WiggleZ sample is difficult, but Table [T] suggests that they 
have similar rest-frame r-band absolute magnitudes. Given 
that typically L , tends to get brighter with redshift (e.g., 
I Wolf et al.ll2003h . this suggests that the WiggleZ galaxies 
are in fact fainter relative to L « in r -band than the SDSS 
sample. However, IWvder et al.l (120071 ) show that the Wig- 
gleZ galaxies occupy the very blue edge of the blue cloud, 
so the fact that both samples are blue does not necessarily 
imply comparable similar formation and evolution scenar- 
ios. With this caveat, if we naively combine the two results, 
we can show that the GI correlations are not significant for 
blue galaxies over a large range of redshifts; we will quantify 
this statement and its cosmological implications in Sec. [6] 
Even if this combination of the two samples is not valid, the 
WiggleZ result is highly useful because of its proximity to 
redshifts used for cosmic shear studies. 

The constraints with the 2SLAQ sample, originally pre- 
sented and interpreted in iHirata et al.l (|2007l ). give a slight 
suggestion of nonzero GI alignments at the 2a level. (For 
reference, we also show the strong positive detection for one 
of several SDSS LRG luminosity bins, which has a higher 
mean luminosity and lower mean redshift.) In contrast with 
the results for red galaxies, the WiggleZ galaxies at the same 
redshift have no such suggestion of nonzero signal, and the 
constraints are significantly tighter because the sample size 
is larger. 

Finally, we presen t a comparison against the results of 
iHevmans et all (|2006d l. who use iV-body simulations to in- 
vestigate the intrinsic alignment signals resulting from differ- 
ent methods of populating the dark matter halos with galax- 
ies (including different ways of aligning the galaxy shapes 
with the dark matter halo or angular momentum vector). 
Our null result is consistent with the results for their "spi- 
ral" model, which is that of a thick disk randomly misaligned 
with the dark matter halo angular momentum vector. The 
mean misalignment angle in that model is 20°. 



6 INTERPRETATION AND COSMOLOGICAL 
IMPLICATIONS 

We now fit two different models to our measured intrinsic 
alignment correlation functions, and investigate the impli- 
cations of our results in terms of the bias in cosmologi- 
cal measurements of the amplitude of matter fluctuations 
from cosmic shear. In order to avoid overly optimistic con- 
straints on model parameters due to inversion of the noisy 
jackknife covariance matrices, we apply the correction de- 
scribed by Hartlap et al. (2007) (equation 17 in that pa- 
per) after inversion. This correction compensates for the fact 
that the inverse of a noisy covariance matrix is not an un- 
biased estimate of the inverse covariance matrix, and corre- 
sponds to multiplication of the inverse covariance matrix by 
(47 — iV r )/47 for 49 jackknife regions and N r radial bins used 
for the fit. Thus, for fits that only use a subset of the radial 
bins, we restrict to the appropriate subset of the covariance 
matrix, invert, and multiply by the appopriate fraction. 

We first fit a power law in transverse separation 
to each of the measured sig nals, in a similar way to 
iMandelbaum et al.1 (|2006bh and IHirata et al] <|2007h . These 
fits have the advantage of being very simple, and allow di- 
rect comparison with previously-published results that use 
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Figure 8. Constraints on the power law amplitude A. Upper 
and lower panels show results for the projected density-shape 
(ws+) and shape-shape (u>^-A correlation functions, respectively. 
The dot-dashed line uses the SDSS Main Blue L4 data from 
M06 and H07; the solid line is for the full WiggleZ sample with 
Ilmax = 60/i _1 Mpc; the dashed line shows results from the 
n m ax = 20h~ 1 Mpc correlation function. In all of those cases, 
the power law slope has been fixed at a = —0.88. The dot- 
ted line shows constraints from our default WiggleZ sample after 
marginalisation over a. 
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Figure 9. Constraints on the power law amplitude A as a func- 
tion of redshift. We use our default line-of-sight range (n max = 
60h~ 1 Mpc), and fix the power law slope at a = —0.88. The up- 
per panel shows constraints from mj s + and the lower panel uses 
From left to right, the points show constraints from SDSS 
Main Blue L4 (cross); WiggleZ z < 0.52 (triangle); and Wig- 
gleZ z > 0.52 (triangle). The horizontal lines indicate the redshift 
range of the observations. We indicate the constraint from the full 
WiggleZ sample (both redshift ranges) with a shaded rectangle. 



the same fitting method; however, they are lacking in physi- 
cal motivation. In order to give fits with physical motivation, 
we fit the unknown a mplitude in a linear a lignment model 
(|Catelan et al.1 l200ll ; IHirata fc Seliakl 120041 ). Using simple 
assumptions, these fits allow us to compare constraints from 
the observed correlation functions and propagate the con- 
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Table 3. 95 per cent confidence limits for power law fits with two 
procedures: varying the power law amplitude A with fixed slope 
a = —0.88, and varying both A and a and marginalising. The best 
fit point in the one or two dimensional parameter space is shown, 
with limits calculated in the one dimensional marginalised space. 
When there are large degeneracies in two dimensions the peak in 
the two dimensional parameter space does not coincide well with 
the peak of the marginalised probability distribution, and can 
even lie outside the one dimensional iso-probability confidence 
limits. For example this occurs in the last line of each table section 
for the parameter a. 



Data 


A 


a 


™s+ 



SDSS Main Blue L4 
WiggleZ n max = 20, all z 
WiggleZ n max = 60, all z 
WiggleZ Il max = 60, 2 < 0.52 
WiggleZ Il max = 60, z > 0.52 
SDSS Main Blue L4 
WiggleZ II max = 20, all z 
WiggleZ n 

max 

= 60, all z 
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0.0000 
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-0.53 
-0.56 



+3.03 
2.62 

+ 3.3(1 
3.43 

+3.74 
3.43 



W++ 



SDSS Main Blue L4 
WiggleZ n max = 20, all z 
WiggleZ n max = 60, all z 
WiggleZ Il max = 60, z < 0.52 
WiggleZ Il max = 60, z > 0.52 
SDSS Main Blue L4 
WiggleZ n max = 20, all z 
WiggleZ II max = 60, all z 



0000+ 0008 
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0040+ 0099 
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straints through to biases on the amplitude of matter clus- 
tering. 

To interpret the observed GI correlation function w g + 
in terms of the correlation between the intrinsic shape and 
the density field (ws+), we assume a linear bias model and 
estimate the bias of the galaxies used as tracers of the den- 
sity field. Values are shown in Table [T] for each of the three 
WiggleZ samples we consider (full s ample and two red shift 
slices), following the methodology of lBlake et al.l (|2009f ). We 
then assume 



= bgWs+ 



(12) 



throughout. A minimum separation of r p > 5/i _1 Mpc is 
used for fits to w g + to avoid non-linear biasing on small 
scales affecting the interpretation of w g + in terms of ws+- We 
use all calculated data points, with no minimum separation, 
when fitting to the ellipticity-ellipticity correlation function 
W++. 



6.1 Power-Law Fits 

We fit a power law 



W5+ = As+ 
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(1 - /bad)' 



(13) 



to the projected density-shape correlation function, and sim- 
ilarly for the II (shape-shape) correlation function w++ . This 
fit is performed separately for each sample under considera- 
tion. Here /b a d is the fraction of objects with bad redshifts, 
as given in Table [T] We do not include a dependence on 
redshift in this subsection, because it would not be mean- 
ingfully constrained by any one dataset alone. We instead 
defer questions of redshift evolution to the following subsec- 
tion. 

We compute the likelihood as a function of A and a on 
a grid using Pr(D|A, a) oc exp(— X 2 /2), using the full covari- 
ance matrix of the data in the \ 2 calculation. We used suf- 
ficiently wide ranges of A and a that the constraints shown 
are not affected. We calculate the resulting amplitudes A 
for two cases: fixed a — —0.88, or marginalised over a after 
allowing it to vary within a wide range. This fixed value for 
a is motivated by line 5 of table 6 of H07, which gives results 
for SDSS LRGs using a minimum separation of 7.5/i." 1 Mpc. 
95 per cent confidence limits were calculated by finding the 
iso-probability level containing 95 per cent of the probabil- 
ity. 

The solid and dashed lines in Fig.[8]show constraints on 
the power law amplitude for the two line-of-sight integration 
ranges considered in this paper, II max = 20 and 60/i _1 Mpc. 
As expected, the larger line-of-sight range results in a weaker 
constraint due to dilution of the signal by noise from uncor- 
rected pairs. The results from the SDSS Main Blue L4 sam- 
ple (M06, H07) are shown for comparison (dot-dashed line). 
The SDSS results give a much tighter constraint on the am- 
plitude of the ellipticity-ellipticity power law than the Wig- 
gleZ samples, and appear close to a delta function at A = 
on the lower panel of Fig. [5] For the ellipticity-density cor- 
relation function, the constraints from SDSS and WiggleZ 
are more similar. We also show the results after marginal- 
ising over the spatial power law coefficient a. These results 
are peaked at A = because of the large range of a values 
allowed if A = 0. 

Fig. shows the 95 per cent confidence ranges as a 
function of sample redshift for the SDSS data and WiggleZ 
redshift subsamples. Since all points are consistent with zero, 
there is no sign of a trend with redshift. However, a strong 
redshift evolution could be ruled out. 

These results are summarised in Table [3] (95 per cent 
confidence limits). We also show results when both the am- 
plitude and scale dependence of the power law are varied. 
Since the amplitude is nearly consistent with zero, the con- 
straints on the power law slope are relatively weak and easily 
consistent with the fiducial value of a — —0.88 taken from 
the SDSS LRG sample. We interpret the w g + constraints as 
intrinsic alignments constraints, and the w++ constraints as 
non-detections of both intrinsic alignments and significant 
shear systematics. 

6.2 Linear Alignment Model Fits 

We now fit a simple but physically motivated intrinsic 
alignment model to the observed correlation functions. We 
us e a variant of the linear alignment model described 
in lCatelan et al.l (|200lh . As originally discussed in Section [5] 
our non-detection of GI correlations constrains the intrinsic 
alignments in the context of this model so that the allowed II 
signals are well below the size of our errors on w++. Thus, 
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for this section we will only use u> g + for constraining the 
linear alignment model. 

The linear alignment model assumes galaxies are 
stretched by an amount proportional to the local curva- 
ture of the smooth ed gravitational p otent ial, and was de- 
veloped further by iHirata fc Selia k (2004) to predict the 
contributions to the lensing power spectraQ This model is 
not expected to be an accurate description o f alignments on 
non-linear scales, and iBridle fc Kind (|2007T l were inspired 
by H07 to insert the non-linear matter power spectrum into 
the model in place of the linear matter power spectrum. 
We refer to this model as the non-linear power spectrum 
linear alignment (NLA) model hereafter. A physically mo- 
tivated mode l for smaller scales based on the halo model 
presented by ISchneider fc Bridle! (|2010T ) gives qualitatively 
similar results to the NLA. Note that while we use the non- 
linear matter power spectrum to describe the density field 
on small scales, we still need the linear bias assumption to 
relate w g + to u>s+ , and therefore cannot use very small scales 
when interpreting this observation (unlike for w++). 

We use the NLA for the rest of this Section, ignoring a ll 
but the first term in equation (16) of lHirata fc Seliakl l|2004h . 
It has a single free parameter, the amplitude Ci, which is 
the constant of proportionality between the galaxy ellipticity 
and the local potential curvature. The model has an inbuilt 
motivated variation with scale and redshift. However, it is 
possible that the constant of proportionality C\ may addi- 
tionally vary with environment and galaxy type, as well as 
redshift. In this paper we therefore allow the amplitude of 
these alignments to vary with an additional free power law 
in redshift, with index pother, to include any variation with 
redshift due to other physics not included in the NLA model. 
The default value is pother = 0, so that the inbuilt redshift 
dependence of the NLA is recovered. Our physical motiva- 
tion for including this additional redshift evolution factor is 
to allow for more complicated effects in galaxy evolution, 
such as mergers and interactions. For example, a model in 
which galaxies align with the local tidal field at early times, 
but gradually decorrelate due to interactions and mergers, 
would be described with a positive ?7othcr that allows for 
significant alignment at early times that is not present at 
the current time . These models are part icularly of interest 
since, as noted in iDrinkwater et all |2010l ), a significant frac- 
tion of the WiggleZ galaxies appear to be interacting, merg- 
ing, or have recently undergone a merger. The full, redshift - 



2 It has been found (Hirata 2010 in prep, an d Joachimi et al 
2010 i n prep.) that the linear alignment model in lHirata fc Seliakl 
has an error in the derivation of its redshift evolution, and 
should be multiplied by a factor of 1/(1 + z) 2 . Thus, the corrected 
version of the linear alignment model corresponds to our zNLA 
model with additional redshift evolution ri ot ^ eI = —2, which as 
we will show in the right panel of Fig. 1121 corresponds to the 
peak of the likelihood of Tfothei (when using the combination of 
low-redshift SDSS data plus WiggleZ). Given that (a) this issue 
was discovered after these calculations were completed, (b) our 
constraints on r] ot ^ er are not very tight, and (c) the constraints on 
contamination to cosmology do not change since we integrate over 
a large range of Mother > we have opted to leave all calculations and 
figures in terms of the lHirata fc S eliak (20o3) formulation of the 
linear alignment model. However, the reader should keep this in 
mind when considering our results in the context of implications 
for linear alignments. 



Table 4. 95 per cent confidence limits and 95 per cent one-tailed 
upper or lower limits for zNLA fits to the GI correlation function 
ui 9 + using two procedures: varying the model amplitude Ci with 
no extra redshift evolution (?7 th er = 0) beyond that already in the 
NLA, or varying both C\ and the extra free power law in redshift 
Mother- d is in units of 5 X lO- 14 (ft 2 Af Mpc- 3 ) -1 The WiggleZ 
results all use n max = 60ft — 1 Mpc. Unless otherwise stated, all 
WiggleZ results use the two WiggleZ redshift bins z < 0.52 and 
2 > 0.52 instead of the "all z" single redshift bin. 



Data 

SDSS Main Blue L4 
WiggleZ all z 
WiggleZ 

SDSS Main Blue L4 + WiggleZ all z 
SDSS Main Blue L4 + WiggleZ 
SDSS Main Blue L4 
WiggleZ all z 
WiggleZ 

SDSS Main Blue L4 + WiggleZ all z 
SDSS Main Blue L4 + WiggleZ 



Ci 

0.84±i;p 

24+ 1 ' 27 

15+ 103 
u - 10 -1.07 

4fi+ 100 

ft oy+0.85 

U-' 3 ' -0.89 

.23 



'/othc 



> 

0.02 
0.02 
0.02 
0.02 



+6.86 
4.58 
+5.26 
3.48 
+ 1.15 
0.43 
+0.99 
-0.41 



> 6.4 

-19.2 
-19.2 
-19.2 
-19.3 

1 n+ 12 ' 6 

1 - u -12.2 

+ 14.3 
12.3 



-19.6 
-17.1 



-1.7^ 



dependent model with the power-law evolution in addition 
to the NLA model is called the "zNLA model." 

Thus, our model for the E-mode power spectrum of the 
density- weighted intrinsic shear, P?j E (Eq.[7j), including this 



additional term is 



-2 

D 2 (z) 



Ps(k) 



1 + 2 
1 + 2piv 



2r)othc 



(14) 



Likewise, the cross-power spectrum between the intrinsic 
shear and density field, Pgjzi (Eq. [4]), is 



(5,-y 



(*) = 



dp 



Ps{k) 



1 + Zpiv 



(15) 



Here, D(z) oc (1 + z)D(z) is the rescaled growth factor nor- 
malised to unity during matter dominati on. For Ps(k), the 
nonlin ear matter power spectrum, we use lPeacock fc Doddsl 
(|l996f ). We use a pivot redshift of z p i v = 0.3 for the redshift 
power law factor. 

We calculate the correlation functions from the pre- 
dicted _rjower_srjectra using Eq. 23 of H07 for w 3 + and Eq. 
10 of Br idle fc Kind (|2007l) . which assumes the ellipticity- 
ellipticity correlation function is equal to its 45 degree ro- 
tated counterpart, i.e. w++ — w X x- These theoretical pre- 
dictions assume that all of the correlation function signal is 
integrated along the line of sight. For the n max = 60ft -1 Mpc 
results, this is a good approximation to the observational 
calculation, but for the n max = 20ft -1 Mpc results this is 
not expected to be the case. The exact amount of signal lost 
by cutting at this shorter line-of-sight integration range will 
depend on modelling of the redshift-space distortions, which 
is beyond the scope of this paper. Therefore for the remain- 
der of this paper we use the more conservative line-of-sight 
integration range n max = 60ft -1 Mpc. 

Before comparing with the data, we average the pre- 
dicted model correlation function over the redshift range of 
the data, to obtain the redshift-averaged theory prediction 
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Figure 10. Constraints on the zNLA amplitude parameter C\. 
The dot-dashed line uses Main Blue L4 SDSS data only; the 
dashed line is from a joint analysis of the two WiggleZ rcdshift 
bins (WiggleZ z < 0.52 and WiggleZ z > 0.52) for n max = GOti" 1 
Mpc; the solid line combines the Main Blue L4 SDSS constraints 
with those from WiggleZ. The dotted line marginalises over addi- 
tional rcdshift evolution beyond that already in the NLA ( r? thcr) 
whereas the other lines all assume no additional rcdshift evolution 
(bother = 0) . 



(zNLAh 



<< NLA) >, - 



dzW(z)w^f LA \z) 



(16) 



Here, the appropriate weight function W(z) is proportional 
to the squared sample redshift distribution, p 2 (z), after di- 
viding out the comoving volume-redshift relation dV c /d;z: 
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Figure 11. Constraints on the zNLA amplitude C\ and an extra 
power law variation with redshift, with index Mother- The up- 
per left and lower right panels show the likelihoods for the two 
parameters separately, where in each case we marginalise over 
the other parameter. The lines on the upper right two dimen- 
sional contour plot contain 68 per cent and 95 per cent of the 
probability. The dot-dashed lines use SDSS Main Blue L4 data 
only; the dashed lines use both the WiggleZ rcdshift bins; and 
the solid lines show the joint constraints from SDSS Main Blue 
L4 and both WiggleZ bins. For comparison we also show results 
using the single WiggleZ dataset for the whole redshift range, 
combined with SDSS Main Blue L4 (dotted line). We use priors 
|Ci/(5 x 10- 14 (h 2 M s Mpc- 3 )- 1 )| < 13 and |ry ot horl < 30 which 
do affect the positions of the dot-dashed and dashed lines for the 
two surveys taken separately. This prior-dependence comes from 
the fact that the SDSS does not cover a wide enough redshift 
range to strongly constrain Mother on its own. We zoom in on the 
range |Ci/(5 X 10- 14 (/i 2 M Q Mpc- 3 )- 1 )| < 2 for clarity. 



for comoving volume V c out to redshift z. The derivation of 
this weight function is given in Appendix [X] essentially it 
comes from the fact that the number of pairs is determined 
by the comoving number density of galaxies, but we are 
integrating over redshift rather than volume. The redshift 
distribution for WiggleZ is taken to be the double Gaussian 
described in Section |3T] When fitting this model prediction 
to the data, we correct statistically for the fraction of bad 
redshifts in the same manner as for the power-law fits: 

^+ = (^ Z + NLA) ), (l-/ bad ) 2 . (18) 

Table [4] shows constraints from each of the two cor- 
relation functions w g + and w++ separately, as well as 
joint constraints obtained by multiplying likelihoods from 
each. Joint constraints are possible now that the under- 
lying model has a single free parameter that propagates 
into each of the two correlation functions. When both 
the amplitude C\ and additional redshift variation power 
law index pother are varied, the constraints depend on the 
prior ranges used for both parameters. We use wide ranges 
|Ci/(5 x lO-^^MoMpcr 3 )- 1 )! < 13, |r?other| < 30, so 
that no constraints are affected by these exact values, except 



those that will never converge no matter how wide the ranges 
are. For the SDSS data alone, since its redshift distribution is 
very narrow and the pivot redshift is well above the redshift 
of the galaxies, the additional redshift power law index is 
unbounded from above, and we show one-tailed 95 per cent 
limits, given the priors. These SDSS constraints on Ci and 
pother are not very meaningful, because the pivot redshift 
was chosen to optimise constraints that use the WiggleZ 
data. However, while the results in terms of Ci and pother 
would appear more meaningful if the calculations used a 
pivot redshift z p i v = 0.1 for the SDSS results, the estimates 
of the bias on erg from cosmic shear surveys are insensitive to 
the choice of pivot redshift. So for consistency we simply use 
the same pivot redshift for all (Ci,»jothor) constraints, and 
emphasise that the numerical values for the SDSS sample 
alone in Tableware not very meaningful. 

Note that u>++ is proportional to the square of the NLA 
amplitude parameter, so the likelihood function is symmet- 
ric about Ci = from these data alone. We do not have a 
strong detection of a non-zero amplitude for any combina- 
tion of w++ data considered, so in the table we show one- 
tailed 95 per cent upper limits calculated from the region 
with Ci > 0. The constraints are much tighter from w g + 
than from w++, which becomes clear when comparing the 
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Figure 12. Left three panels: Maximum and minimum bias on ag allowed within the 95 per cent confidence limits of the fits shown 
in Fig. 1111 at each value of the additional rcdshift dependence parameter, rather- Line styles are identical to those in Fig. 1111 if a line 
stops at some maximum value of 7? ot hen it means that all higher values of »7 ot her are ruled out by the data at this confidence level. We 
assume a CFHTLS-like survey and carry out a Fisher matrix analysis assuming all cosmological parameters are known except erg. The 
left hand panel artificially sets the II contribution to the cosmic shear signal to zero to show the impact of Gl contamination alone; the 
second panel shows the impact of II contamination alone; and the third panel shows the expected full intrinsic alignment contamination, 
including both Gl and II contributions. Right panel: From Fig. 1111 the likelihood as a function of Mother i f° r easy comparison with the 
other panels. This plot can be used to judge how likely any of the erg biases in the other panels are. 



joint constraints (bottom section of the Table) with the con- 
straints from w g + alone (top section of the Table). For the 
remainder of this section we show only the joint constraints. 
Note that this finding is perfectly consistent with the power 
law results shown in the previous subsection, for which the 
power law amplitude was more tightly constrained from the 
w++ data than the w g + data. It arises because the NLA 
prediction for w++ is significantly smaller than the w++ er- 
ror bars. From the Table we see that the results for both 
the amplitude and the additional power law index are all 
consistent with zero. 

Fig. [10] shows constraints on the NLA amplitude pa- 
rameter Ci assuming the standard NLA redshift evolution 
(pother = 0). SDSS Main Blue L4 (dot-dashed line) has some- 
what less constraining power than WiggleZ (dashed line). 
This result may at first seem in contradiction with the power 
law constraints, for which SDSS gave a much tighter con- 
straint on the power law amplitude than WiggleZ. However, 
it makes sense because the constraints are consistent with 
Ci = and the NLA predictions for SDSS have absolute 
values closer to zero than for WiggleZ. They are therefore a 
similar size relative to the error bars and thus give similar 
constraints on C\. (See also a similar effect, discussed above, 
regarding the relative constraining power of w g + and w++.) 
The joint constraint (solid line) is consistent with zero. We 
also show the constraint on the amplitude parameter after 
marginalising over the additional redshift dependence (dot- 
ted line). This procedure gives a sharp spike at C\ = due 
to the large range in rather allowed at zero amplitude. 

Joint constraints on the amplitude and additional red- 
shift power law index are shown in Fig. 1111 The SDSS data 
alone can only place a lower limit on the additional redshift 
dependence (dot-dashed line) . Therefore the limits shown in 
the Table depend completely on the prior ranges used for 
Ci and pother- For the particular prior ranges used here, the 
lower limit is suggestive of strong redshift dependence (e.g. 
pother > 7.7 for the joint w g+ + w ++ dataset). This is to be 
expected from the green dot-dashed lines in the lower right 



panel of Fig. 1111 but should be taken with a grain of salt 
due to the dependence on priors. 

The WiggleZ data alone (dashed line) place both an 
upper and weak lower limit on the redshift variation due 
to the redshift range spanned by the data. The combined 
result (solid line) constrains the power law index to lie in on 
a relatively narrow region, but this region still has a width 
of around 7 powers in redshift. The dotted line shows for 
comparison the joint result using only the combined WiggleZ 
sample, which is similar but slightly weaker than the joint 
result using both subsamples, as expected. 

6.3 Cosmological interpretation 

In this subsection, we estimate the bias in the measured 
linear theory present day amplitude of density fluctuations 
as, if intrinsic alignments due to blue galaxies were ignored 
in a cosmic shear analysis. While we expect this bias to 
be consistent with zero, since we find w g + and w++ to be 
consistent with zero, we would like to find how tightly we 
can constrain the bias. For this purpose, a Fisher matrix is 
calculated using a single survey redshift b in (no tomogra- 
phy) w ith the redshift distribution given in iBeni amin et al.l 
(|2007h for the Canada- France-Hawaii Telescope Legacy Sur- 
vey (CFHTLSfl Th is is propagated into the bias on param- 
eters using eq. 21 of lHuterer et al.l (|2006t ). where in practise 
we only consider ag and assume all other parameters are 
fixed. Note that the sky area and number density of galax- 
ies drop out of this calculation. 

We assume for simplicity that the galaxies used for cos- 
mic shear are the same as both the WiggleZ and L4 galaxies, 
though as previously discussed these two samples are not 
necessarily comparable in formation history. Consequently, 
we will also consider how much the bias can be constrained 
using WiggleZ and L4 blue galaxies separately. In practise, 
the CFHTLS or other comparable cosmic shear surveys will 

3 http : //www . cf ht .hawaii . edu/Scienc e/CFHLS/ 
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incl ude some red gal axies (roughly 20 per cent of the sam- 
ple, [WolfeF^aO [200^) which likely have a stronger intrinsic 
alignment signal (MOB, H07), and will also include fainter 
blue galaxies which may have a weaker signal. These flux- 
limited cosmic shear samples will tend to be dominated by 
L* galaxies such as those in the SDSS L4 sample and in the 
WiggleZ survey (which spans a range of luminosities but 
has a mean around L, at the mean redshift). In terms of 
the colour distribution, the SDSS Main L4 blue sample con- 
tains galaxies with colours spanning the entire blue cloud, 
whereas the WiggleZ survey contains the bluest ~ 10 per 
cent of the blue cloud. Consideration of more complicated 
modelling of the intrinsic alignment amplitude as a func- 
tion of luminosity, colour and redshift is beyond the scope 
of this paper. However, the impact of red galaxy intrinsic 
alignments on such a survey was already estimated by H07 
using the measured signals from SDSS and 2SLAQ. 

To produce Fig. 1121 we consider one value of the addi- 
tional redshift power law index pother at a time, for all pother 
values allowed within the 95 per cent confidence region of 
Fig. 1111 We then calculate the bias in erg for each value of Ci 
allowed within the 95 per cent confidence region of Fig. [11] 
for that value of pother , and find the maximum and minimum 
bias value. We plot these bias values as a function of pother 
in Fig. 1121 We repeat this procedure for each of the dataset 
combinations shown in Fig. 1111 

We consider separately the case where the II contri- 
bution to the cosmic shear power spectra is zero and the 
only contamination comes from the GI term (left panel of 
Fig. I12[> . Similarly we consider the II-only case in the middle 
panel of Fig. [T2j While such configurations are not possible 
within the NLA, since Ci is the same parameter figuring 
into both the GI and II correlations, this separation gives us 
some additional physical understanding of the constraints. 
Finally, the right panel of Fig. H2l shows the total bias on erg 
taking into account both contributions. 

Starting in the left panel of Fig. 1121 we see that we 
cannot rule out large positive or negative GI contamina- 
tion when using SDSS alone, because (as we have already 
seen) its short redshift baseline makes it impossible to rule 
out significant positive pother- WiggleZ is able to place more 
stringent constraints due to its higher mean redshift and 
much broader width. The combination of the two surveys is 
able to effectively narrow the constraints on the GI contam- 
ination both for positive and negative pother- In the second 
panel, we see similar effects due to the survey redshift dis- 
tributions. As expected for II contamination, the bias on 
the amplitude of fluctuations is always positive. In the third 
panel, we show the combined intrinsic alignment effects. The 
right-most panel, which appeared in Fig. 111! can be used to 
evaluate the likelihood of the erg biases in the other panels. 
For example, since the combined SDSS and WiggleZ sam- 
ples (solid line) constrain redshift evolution on top of that 
predicted by the linear alignment model, we can see that 
the large as bias for pother > 5 is quite unlikely. However, 
for the SDSS alone (dot-dashed line), large positive 7/ ther 
(corresponding to significant alignments at higher redshift) 
cannot be ruled out, and thus the as biases at large pother 
also cannot be ruled out in the absence of external reasons 
to discount large father- 

If we believe the NLA is the sole (or main) contribu- 
tor to the blue galaxy intrinsic alignments with no extra 



redshift evolution (7/ ot her = 0), then the range of bias Acrg 
allowed within the 95 per cent limits on C\ and rather are 
are -0.03 < Aa 8 < 0.02 from SDSS, -0.03 < Acrg < 0.03 
from WiggleZ, and -0.03 < Acrg < 0.02 from the two sur- 
veys combined. When using both surveys together, the con- 
straints are similarly powerful for models with increasing 
alignments at lower or higher redshift (i.e. varying pother)- 
We have already ruled out models with pother > 7, and can 
rule out large biases for all allowed values of pother- When 
using SDSS alone, we cannot rule out very large Acrg > 0.1 
due to blue galaxies if rather > 3, which could be problem- 
atic if there is some physics that causes such strong red- 
shift evolution. When using WiggleZ alone, our constraints 
on the erg bias in cosmic shear measurements dominated by 
blue galaxies are roughly ±0.04-0.06 (95 per cent CL) for 
Mother < 0, which is weaker than when we include SDSS blue 
L4 galaxies, but still sufficient to rule out intrinsic alignment 
systematics that are comparable to the statistical error for 
current cosmic shear surveys. 

The main caveats regarding these constraints relate to 
the nature of the samples used and the models used to in- 
terpret the data. In particular, as already mentioned, the 
WiggleZ and SDSS blue L4 galaxies may have different for- 
mation histories, in which case the comparison of their re- 
sults may not be meaningful. Furthermore, we neglect the 
red galaxies, for which constraints have already been placed 
using a different procedure in H07. Finally, we have not at- 
tempted to interpret the data in light of other models for 
intrinsic alignments, of which several exist in the literature; 
use of other models, or additional redshift dependence that 
is poorly modelled by a power-law, may change the pro- 
jected intrinsic alignment contamination from the numbers 
shown here. 



7 CONCLUSIONS 

In this paper, we have placed the first direct observational 
constraints on the intrinsic alignments of blue galaxies at 
intermediate redshift (z ~ 0.6), using the WiggleZ spectro- 
scopic redshifts and galaxy shape measurements from SDSS. 
We followed a comparable procedure as has been used be- 
fore in SDSS at low redshifts (z ~ 0.1- 0.3) for blue and 
red ga laxies in iMandelbaum et all i|2006ri ) and lHirata et all 
l|2007h . This procedure relies on finding pairs of galaxies that 
are physically associated in terms of their three-dimensional 
separation, and calculating the correlation between their 
shapes, and between the shape of each galaxy with the line 
connecting their positions on the sky. 

Our result was a null measurement for the full Wig- 
gleZ sample and for two redshift subsamples. This null mea- 
surement can in turn be used to constrain parameters of 
physically-motivated intrinsic alignment models, and to con- 
strain the contamination of cosmic shear observations due to 
intrinsic alignments of galaxies that are comparable to this 
sample. We have found that if we assume a model involving 
linear alignment with the smoothed local density field, then 
we can constrain the intrinsic alignment contamination for 
a CFHTLS-like survey dominated by WiggleZ-like galaxies 
to be small enough that erg is biased by an amount that 
is smaller than the statistical errors. If we allow additional 
power-law redshift evolution in these alignments on top of 
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the redshift evolution that is encoded in the linear alignment 
model, then we see that the constraints for pother | < 2-4 do 
not significantly weaken, and the models with very large 
pother can be ruled out because the data cover a fairly long 
redshift baseline (but see footnote 2 in section ROl for a cau- 
tionary note on interpreting these pother values). Combina- 
tion with low-redshift SDSS results, which may be valid if 
the UV-selected WiggieZ galaxies have comparable forma- 
tion histories to optically-selected L* blue cloud galaxies in 
SDSS, allows for tightening of these constraints, particularly 
due to the ability to rule out models with strongly increasing 
intrinsic alignments at low redshift. 

As previously noted, theoretical models of intrinsic 
alignments are currently poorly constrained by the data, 
so direct measurements of these alignments are necessary to 
estimate how serious the alignments are for current and fu- 
ture cosmic shear surveys. These observations of blue galax- 
ies at intermediate redshift fill in an important gap in our 
knowledge of intrinsic alignments. Our constraints were pre- 
dominantly phrased in terms of the linear alignment model. 
However, we note that our null measurement if II alignments 
could also be used to constrain other intrinsic alignments 
models, such as the quadratic alignment model that predicts 
no GI alignments but potentially significant II. 

While blue galaxies tend to dominate cosmic shear sam- 
ples, strong alignments of red galaxies at intermediate-high 
redshift could still be a significant contaminant. As a re- 
sult, it will be important to obtain similar constraints of 
intrinsic alignments of red galaxies, which are poorly con- 
strained for z > 0.4. Future work with our measurements of 
intrinsic alignments in WiggieZ could also focus on consider- 
ing other physically-motivated intrinsic alignments models, 
and on combining pre-existing constraints for red galaxies at 
low redshift with our constraints for blue galaxies to come 
up with an estimate of intrinsic alignment contamination 
to cosmic shear surveys with a more realistic blue plus red 
galaxy sample. 

Several methods have been proposed to remove the in- 
trins ic alignment signal from future cosmic shear surveys 
(e.g. iKingj 20051; iBridle fc Kindl2007l; Ijoachimi fc Schneider! 
20081 . l2009l ; IZhang!l2008l ; lBernsteinll2009l ; |joachimi fc Bridie] 
20091 ). In general, these methods rely on using some of the 
weak lensing signals (auto- and cross-correlations) to con- 
strain parameters of the intrinsic alignment models, result- 
ing in a loss of cosmological information. Our measurements 
in this paper using the WiggieZ dataset will allow for the 
placement of stronger priors on the intrinsic alignment mod- 
els, and therefore minimise this loss of cosmological infor- 
mation, preserving the cosmological constraining power of 
future datasets. 
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APPENDIX A: WEIGHT FUNCTION FOR 
REDSHIFT- AVERAGING OF THEORETICAL 
SIGNAL 

Here, we derive the weight function to use when averag- 
ing the theoretical signal (such as the NLA model) over the 
broad redshift distribution of the WiggleZ sample, for com- 
parison with the observed signal. For simplicity, the calcu- 
lation is done in terms of a simple estimator for the galaxy 
autocorrelation £ 93 ; however, the conclusions can clearly be 
applied to other pair statistics such as the intrinsic align- 
ment cross-correlation £ 9 +. We phrase all calculations in 
terms of the following estimator: 



DD 
~RR 



(Al) 



This estimator, which uses data-data and random-random 
pairs, is similar to our estimator for £ s +, Eq. JS} except that 
the latter also requires data-random pairs. 

When calculating the correlation function, we find all 
galaxy pairs in a particular bin in (r p , II) with volume Vbin, 
and accumulate them without any regard for their redshift. 
We then make an estimated £(r p ,II) which presumably is 
some weighted function of z, i.e. 



f= / dzt(z)W(z). 



(A2) 



Our goal is to derive the functional form of W(z). 

The key point to understand is that we average over 
all pairs, and thus the weight function is determined by the 
number of random-random and data-data pairs RR and DD 
that belong in this bin in r p and II. The former is sim- 
ply proportional to the comoving volume density of galaxies 
multiplied by the volume factors, and the latter has an ad- 
ditional factor of 1 + £: 



RR <x[$z) VbindK 
DD oc (1 + £)RR 



(A3) 
(A4) 



for a comoving volume V c . 

When we calculate £ in one big redshift bin, ignoring 
the redshift of any real or random pair, we effectively sum 
(integrate) RR and DD over our whole redshift range before 
using our estimator to get £. Thus, 



JdDD 
JdRR 



1 = 



JdV c ^(dN/dV c ] 
J dV c (dN/dV c )- 



(A5) 



We can then transform the integrals so they are over red- 
shift. In addition, we define dN/dz = Np(z) where JV is the 
total number of galaxies, and p(z) is integrates to 1. 
As a result, Eq. (|A5|) becomes 



jp\z)(dz/dV c )i(z)dz 
j p 2 (z)(dz/dV c )dz ' 



(A6) 



The single factor of dz/dV c comes from the squared num- 
ber density of galaxies requiring a (dz/dV c ) 2 , and the dif- 
ferential volume element over which we integrated before, 
dV c , becoming dz(dV c /d,z). Thus, comparison of Eqs. (|A6|I 
and (|A2|I allows us to identify the redshift weight function 



W(z) 



p 2 (z)(dz/dV c ) 
Jp 2 (z)(dz/dV c )dz 



(A7) 



The reason this is not simply p 2 (z) is that the number of 
pairs scales with the volume number density, so when inte- 
grating over redshift we get factors of the rate of change of 
comoving volume with redshift. 
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